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ABSTRACT 

Context. The ground-based campaigns on S Set stars revealed richer and richer frequency spectra of these opacity-driven pulsators, 
thanks to continuous improvements in their exploitation. It has also been suggested that the detection of a wealth of very low amplitude 
modes was only a matter of signal-to-noise ratio. Access to this treasure, impossible from the ground, is one of the scientific aims of 
the space mission CoRoT. 

Aims. This work presents the results obtained on HD 50844, the only 6 Set star observed in the CoRoT Initial Run (57.6 d). The 
aim of these CoRoT observations was to investigate and characterize for the first time the pulsational behaviour of a 6 Set star, when 
observed at a level of precision and with a much better duty cycle than from the ground. 

Methods. The 140,016 datapoints were analysed using independent approaches (SigSpec software and different iterative sine-wave 
fittings) and several checks performed (splitting of the timeseries in different subsets, investigation of the residual light curves and 
spectra). A level of 10~ 5 mag was reached in the amplitude spectra of the CoRoT timeseries. The space monitoring was complemented 
by ground-based high-resolution spectroscopy, which allowed the mode identification of 30 terms. 

Results. The frequency analysis of the CoRoT timeseries revealed hundreds of terms in the frequency range 0-30 d~' . All the cross- 
checks confirmed this new result. The initial guess that S Set stars have a very rich frequency content is confirmed. The spectroscopic 
mode identification gives theoretical support since very high-degree modes (up to ^=14) are identified. We also prove that cancellation 
effects are not sufficient in removing the flux variations associated to these modes at the noise level of the CoRoT measurements. The 
ground-based observations indicate that HD 50844 is an evolved star that is slightly underabundant in heavy elements, located on 
the Terminal Age Main Sequence. Probably due to this unfavourable evolutionary status, no clear regular distribution is observed in 
the frequency set. The predominant term (/i =6.92 d~') has been identified as the fundamental radial mode combining ground-based 
photometric and spectroscopic data. 

Conclusions. The CoRoT scientific programme contains other 8 Set stars, with different evolutionary statuses. The very rich and 
dense frequency spectrum discovered in the light curve of HD 50844 is the starting point for future investigations. 

Key words. Stars: variables: 5 Set - Stars: oscillations - Stars: interiors - Stars: individual: HD 50844 
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1. Introduction 

A huge effort has been made in the past decades to observe 
5 Set stars from the ground. These opacity-driven pulsators 
show many excited modes with amplitudes that can be reached 
by photometry from the ground. After several results obtained 
by di fferent teams performing single-site observations dPorettil 
2000), the most promising targets have been monitored by 
means of multisite ca mpaigns , main ly in the framework of the 
Delta Scuti Network dBregeil |2000|) and of the STEPHI net- 
work dMichel et all I2000I) The most observed 6 Set star from 
the ground is FG Vir: 75 frequencies have been detected with 
amplitudes down to 0. 2 mmag thanks to a multisite monitoring 
spanning several years dBregeret al.U2005l) . 
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The first S Set st ar monitored from s pace was 8 2 Tauri, with 
the WIRE satellite dPoretti et al.L |2002|) . WIRE data, collected 
in only 19 days, increased the number of the known frequencies 
from 5 to 12. More recently, the Canadian Space Agency mission 
MOST detected i n the light curve of the 5 Set star HD 209775 
dMatthewsl 120071) roughly the same number of frequencies as in 
FG Vir. The MOST observations, concentrated in two runs last- 
ing 14 and 44 days and separated by about 1 year, gives a taste of 
the new horizons that can be opened by a space mission, which 
monitors 5 Set stars in a continuous way for several months at a 
time. 

The asteroseismologic scientific case of the CoRoT 
(COnvection, Rotation and planetary Transits) space mission 
gives such an opportunity for long-term monitoring of the 
same f ield dBaglin et all 120061: [Michel et alll2006HCatala et al.L 
120061) . lAuvergne et all d2009b describe the CoRoT in-flight per- 
formances. Some preparatory work has been undertaken to 
identify potentially interesting targets in the selected observing 
zones and, among them, several 5 Scuti stars have been pro- 
posed dPoretti et al.L I2003L 120051) . After the successful launch 
in December 2006 and the in-orbit calibrations, CoRoT ob- 
served its first field in February and March, 2007, i.e., the first 
Initial Run (IR01). One of the new 6 Set stars discovered in the 
preparatory work, HD 50844, was among the 10 asteroseismo- 
logic targets of the IR01. It is a slightly evolved star, close to 
the Terminal Age Main Sequence. The physical parameters de- 
rived from Stromgren photomet ry are My-1 . 31, T e ff=7500 K, 
logg=3.6, and [Fe/H]=-0.4 dex dPoretti et al.L 120051) . In this pa- 
per we review our understanding of 5 Set stars and unveil a new 
"look" of their frequency spectra, thanks to the intensive CoRoT 
monitoring. We also try to answer some open questions on mode 
excitation in the pulsators located in the lower part of the insta- 
bility strip. To give as complete as possible answers, we also ob- 
tained high-resolution spectroscopy and multicolour photometry 
from ground. In particular, the study of line profile variations is 
a complementary approach, and essential for the identification 
of the excited modes. 

2. The CoRoT data 

The IR01 started on February 2 nd , 2007 and finished on March 
31 st , 2007 (AT=57.61 d). The exposure time in the CoRoT as- 
teroseismology channel is 1 sec. For the analysis, we used the 
reduced N2 data rebinned at 32 sec and we only considered the 
140,016 datapoints for which no problem (i.e., flag=0) were re- 
ported. The amplitudes of the background orbital variations are 
kep t at a very smal l level by the great effectiveness of the baf- 
fle dAuvergne et al.L 120091) . The subsequent rejection of the un- 
certain points strongly minimized the orbital effects. The light 
curv e was detrended with a linear fit to remove the effect of age- 
ing dAuvergne et all 120091) . The 32-sec time binning provided 
an oversampling of the light curve coverage. Therefore, to gain 
in CPU time and to reduce the noise level, we grouped the orig- 
inal data into new bins of four consecutive measurements, thus 
obtaining 32,433 datapoints. 

The spectral window is shown in Fig. [1] The rejection of 
the flagged points resulted in a slight enhancement of the am- 
plitude of the orbital frequency, since the bad measurements 
occurred mostly when the satellite crossed the South Atlantic 
Anomaly (SAA). This crossing occurred twice in a sidereal day 
and gave rise to an additional alias peak. In the 0-50 d -1 inter- 
val, there are four structures in the spectral window. The first 
is composed of a peak at 2.006 d -1 (4% of the amplitude) and 
of a lower one at 4.007 d -1 (0.7% of the amplitude). The sec- 
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Fig. 1. Top panel: the spectral window of the CoRoT data on 
HD 50844, observed in the Initial Run of the mission. Middle panel: 
the amplitude spectrum of the original CoRoT data. Bottom panel: the 
power spectrum (in units of the normalized spectral power) of the line 
profile variations detected in the ground-based spectroscopic data. 



ond and most prominent structure is around the orbital frequency 
/ s =13.972 d _1 : f s (3.0% of the amplitude), f s +l d" 1 (7.5%) and 
f s —l d -1 (8.5%). Other structures occur around 2f s (7.0% of the 
amplitude), flanked by lower peaks at2/ s + l d -1 and2/,±2d -1 , 
and around 3/ s (5.0% of the amplitude), flanked by lower peaks 
at3/ s ±l d 1 and 3/,±2d _1 . 

Several app roaches were used to process the CoRoT data, 
e.g., Period04 ( Lenz & Bregerl |2005[). itera t ive si ne-wave fit- 
ting dVanicekL 119711) and SigSpec dReegenl 120071). Moreo ver, 
we refined the classical fitting proposed by I Vanf£ekl d 1 97 lb for 
the analysis of the CoRoT data on classical pulsators with a 
Levenberg-Marquardt algorithm. It uses the frequency domain 
as a first guess for the fi tting of sinusoids in the time domain 
dLefevre & Chendl2009h . 

Although small differences in the values of the detected fre- 
quencies have been found (see Sect. [3]), the main results de- 
scribed here are independent from the approach used. The pres- 
ence of the predominant term suggested in the preparatory work, 
based only on gro und-based data (/i=6.92 d -1 , see Fig. 5 in 
iPoretti et al.L 120051) . has been immediately confirmed. The alias 
structures of this term (e.g., f\ + / s =20.89, f\ + 2/ s =34.86, and 
f\ + 3/s =48. 8 3 d _1 ) are visible in the amplitude spectrum of the 
original data (Fig. Q] middle panel). Both methods were very 
effective in removing the alias structures together with the real 
peak; actually, no relevant term is found at the alias frequencies 
of f\ . Taking the excellent quality of the N2 data into account, 
there is little influence of the satellite's orbital period on the de- 
tected frequencies; in particular, those specifically discussed in 
the paper are free of any aliasing or orbital effect. 
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Fig. 2. Examples of the light curve of HD 50844 observed by CoRoT. The lower panel is a zoom-in of a region of the upper panel. 



3. The frequency content of the CoRoT data 

To have a reference frame for interpreting of the results for 
HD 50844 (A2, 7=9.09), we also considered the data of 
HD 292790 (F8, V=9AS), which was observed by CoRoT in 
the same IR01. Our frequency analysis shows that the luminos- 
ity of HD 292790 is modulated by the rotation in a simple way. 
The amplitude spectrum shows the rotational frequency, its har- 
monics and the satellite frequencies (Fig. [3] top panel). For most 
part of the spectrum, i.e., from 10 to 100 d _1 , the noise level 
is distributed in an uniform way and is very low, namely below 
0.01 mmag. For / < 5 d~\ where the modulation terms and the 
long term drift are concentrated, the noise level increases slightly 
increases. 

In contrast, the amplitude spectrum of HD 50844 appears to 
be very dense (Fig.Q] middle panel, and, in a more striking way, 
Fig. [3] second panel). At first glance, it is clear that the ampli- 
tude spectrum of HD 50844 is only for / > 50 d~' as flat as that 
of HD 292790. We can calculate a noise amplitude of 7.5 x 10~ 6 
mag. It is straightforward to deduce that the signal is concen- 
trated in the / < 30 d~' domain, though there is no predominant 
peak region after identification of 250 frequencies (residual rms 
1.2 mmag). At that point, the average peak height (0.10 mmag) 
for / < 30 d -1 is still 10 times higher than for / > 50 d" 1 . A 
huge number of terms is necessary to reduce the amplitude spec- 
trum of HD 50844 at the expected noise level. The peak height 
gets progressively lower after the identification of 500, 750, and 
1000 terms (third, fourth, and fifth panels in Fig. [3] respectively; 
see also the second column of Table [TJ, but the residual rms 
decreases very slowly (third column of Table [T). Again, many 
low-amplitude peaks are necessary to reduce the residual rms: 
after detection of 500 frequencies, other 500 frequencies with 
amplitudes between 0.06 and 0.03 mmag are needed to reduce 
the residual rms by a factor of 0.60%. 

The distribution of the frequencies offers some clues about 
the nature of the signal. We calculated the distribution of all the 
frequencies, of the first 300 in order of detection and of those 



Table 1. Summary of the frequency detection. 



Detected 
frequencies 


Amplitude Residual rms 
[mmag] 


SNR 


sig 







13.63 






1 


15.45 


8.14 


2060.5 


4526.4 


20 


0.64 


2.85 


85.5 


175.6 


60 


0.28 


2.03 


37.4 


66.3 


100 


0.20 


1.74 


26.7 


46.9 


200 


0.12 


1.34 


16.0 


29.7 


300 


0.09 


1.11 


12.0 


23.4 


500 


0.06 


0.83 


8.0 


18.3 


750 


0.04 


0.61 


5.3 


15.7 


1000 


0.03 


0.48 


4.0 


12.6 



detected after the 300 th rank (Fig. 0). The distribution of all the 
frequencies is almost flat in the 0-20 d _1 region, followed by a 
slow decline until 30 d A This tail is absent in the distribution of 
the first 300 frequencies, and there is a sharp step at 20 d A The 
20-30 d~' region is populated by the frequencies detected after 
the first 300. 

It is worthwhile investigating how the detection of so many 
frequencies modifies the light curves. Figure [5] shows the light 
curve of HD 50844 at different stages of prewhitening. The 
multiperiodic nature of the original data (first panel) is obvi- 
ous. For instance, the sine wave of the predominant frequency 
is almost cancelled by interference with other terms between 
HJD 2454159.5 and 2454159.6. Variability is still evident af- 
ter removing 100 frequencies. This is an important step for- 
ward compared to previous ground-based and space observa- 
tions. Also long-term oscillations are noticed at this level. After 
removing 250 frequencies, the residual light curve shows oscil- 
lations whose periods are around 0.05 d (~ 20 d~'; third panel). 
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Fig. 3. For comparison purposes, the top panel shows the amplitude 
spectrum of the rotational variable HD 292790. This star has been ob- 
served by CoRoT simultaneously to HD 50844. The amplitude spectra 
of HD 50844 shown in the other panels indicate the extreme richness of 
the pulsational spectrum of this S Set star: up to 1000 peaks are needed 
to make the residual spectrum of HD 50844 comparable to the original 
one of HD 292790. 



After removing even more frequencies, the oscillations progres- 
sively become more rapid and periods about 0.03 d (~ 33 d~') 
are visible in the fourth and fifth panels of Fig. [5] This is another 
reflection of the distribution tail observed in the third panel of 
Fig. [4] It is quite evident that oscillations are still present after 
subtraction of 1000 terms. As shown in the last two columns of 
Table Q] after 1000 frequencies the SNR is around 4.0 and the 
Spectral Significance parameter ("sig") of the SigSpec method 
is around 12.6. The former is considered the acceptance limit in 
the frequency search algorithms such as Period04, while the lat- 
ter is still above the threshold of the SigSpec method, i.e., 5.46. 
In a conservative way, we decided to stop the frequency search 
at this point, as it becomes more and more difficult to resolve 
new terms within our frequency resolution. 

We considered it important to check the reliability of the fre- 
quency values determined until here. To do so, we compared 
the results from SigSpec and from the iterative-sine-wave fit- 
ting, considering an error range of 0.005 d in frequency (i.e., 
a width 3.4 times narrower than 1/AT). The cross-checks show 
an almost complete agreement on the values of the first 500 fre- 
quencies: 498 frequencies were found by both methods, whose 
438 (88%) have an amplitude within the 20% range of tolerance. 
In particular, both methods correctly identified the first 462 fre- 
quencies within an interval of few units of 10^ 4 d _1 . These fre- 
quencies are listed in Table 2, which is available in electronic 
form at the CDS. The solutions of the two methods slowly di- 
verge in the 500-1000* range because the differences will ac- 
cumulate in the iterative process of the frequency extraction. 
Notwithstandining, the agreement is still satisfactory when con- 
sidering all the 1000 frequencies: 953 were found by both meth- 
ods, whose 805 with an amplitude within the 20% range. 
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Fig. 4. Distribution of the frequencies. Top panel: histogram of all the 
1000 frequencies. Middle panel: histogram of the first 300 frequencies 
in order of detection. Bottom panel: histogram of the frequencies de- 
tected after the first 300. 



Data - 100 f 



Data - 250 f 



Data - 500 f 



Data - 750 f 



Data - 1000 f 



4159.1 4159.2 4159.3 4159.4 4159.5 4159.6 4159.7 4159.8 4159.9 
HEL J.D.- 2450000 

Fig. 5. The behaviour of the light curve of HD 50844 in the original 
data (top panel) and after subtracting the detected terms at different 
steps, where the time axis spans 1 d. 



We also performed another independent check. We created 
two subsets by considering the first and the second halves of the 
CoRoT timeseries. Again, we found a strong match between the 
frequencies obtained from the two subsets. However, we had to 
take the degradation of the frequency resolution into account and 
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had to relax the error range. Therefore, the same analysis was re- 
peated considering two subsets spanning the same time baseline, 
i.e., 50 days. The first subset was composed of the data in the 
time intervals 0.0-10.0 d and 30-50 d (0.0 d is the time of the 
first measurement) and the second subset consisted of the data 
in the intervals 10-30 d and 50-58 d. The two subsets have 823, 
756, and 501 frequencies (out of 1000) in common considering a 
maximum separation of 0.025 d _1 , 0.020 d -1 and 0.010 d _1 , re- 
spectively. Such large fractions of coincident terms suggest that 
they are more likely intrinsic to the star rather than due to tran- 
sient instrumental effects. It is possible that some of the lower 
peaks originate from changes in amplitude or phase, or from 
other effects intrinsic to the star. In any case, the main conclusion 
of our analysis and related checks is that hundreds of terms - at 
least up to 1000 - are needed to explain the light variability of 
HD 50844. The detection of so many independent terms in the 
light variability of a 6 Set star is a totally new result. 

We also investigated the constancy of the highest amplitude 
terms, i.e., those listed in Table [3] We subdivided the CoRoT 
timeseries into five subsets, each spanning 11.6 d, and we cal- 
culated a least-squares fit of each subset. The amplitude of 
the predominant term f\ term is the most interesting to be in- 
vestigated, since a low percentage of variability will immedi- 
ately result in a detectable effect. The amplitude values span a 
full range of 0.10 mmag, with a formal error on each value of 
+0.05 mmag. Therefore, the amplitude of the /] term can be con- 
sidered constant. No significant amplitude variability has been 
detected in the other terms, except for f&,fi\, and fa. These 
terms show evident amplitude and phase variations. This fact is 
a known effect induced b y the interference of terms clo ser than 
the frequency resolution dBreger & Pamvatnvkhl 120061) . A sec- 
ond investigation of the temporal behaviour of the amplitudes 
was performed usi ng a time-frequency analysis (described in 
iBaudin et al.L 1 1994b . This analysis is based on a Morlet wavelet 
with an adjustable frequency (and thus, time) resolution. If this 
resolution is greater than the frequency difference between two 
neighbouring terms, the power variation shows the beating de- 
scribed above. If the resolution is smaller than the frequency dif- 
ference, the amplitude of the terms can be evaluated. We did not 
found any trace of variability of the amplitudes by performing 
this test on the terms considered above. The constant behaviour 
of the amplitudes is illustrated for the terms f\ and in Fig. [6] 
the small variations observed are compatible with the influence 
of the noise. 

Linear combinations of the terms having the highest am- 
plitudes are commonly found in multiperiodic radial pulsators 
(Cepheids, RR Lyr, high amplitude 6 Set stars), an d are also ob- 
serve d in low amplitude <5 Set stars, e.g., FG Vir dBreger et al.L 
120051) . In the solution of HD 50844 we find some frequencies 
with a high rank of detection that could be explained as a lin- 
ear combination of a frequency with a low rank and the first, 
predominant term f\ (Table[3]l. However, we cannot rule out that 
these high-rank frequencies are actually independent modes, ex- 
cited by resonance mechanisms or even intrinsically excited. 

4. Ground-based observations 

4.1. High-resolution spectroscopy: line profile variations 

The spectroscopic observations were completed from January 2 
to 28, 2007 with the FEROS instrument mounted at the 2.2-m 
ESO/MPI telescope, La Silla, Chile. The spectral resolution of 
FEROS is R ~ 48, 000. We obtained 232 spectra in 14 nights. 
The exposure time was set to 900 sec and the signal-to-noise ra- 




Fig. 6. Time-frequency analysis of the modes 
/i=6.92 d-'=80.1 fiHz (top) and / 4 =12.85 d-' = 148.7 juHz 
(bottom) with a frequency resolution of 2 and 1 //Hz, respec- 
tively. 



tios (SNRs) ranged from 150 to 210. The spectra were reduced 
using an improved version of the s tandard FERO S pipeline, writ- 
ten in MIDAS and developed bv lRainerl d2003h . In addition to 
the individual line profiles, we considered the mean profiles 
computed using the least-squares deconvolution (LSD) method 
dDonati et al.L 1 1997b . When calculating the LSD profiles, the 
spectral region between 4150 and 5800 A was used, taking care 
to omit the intervals containing the Hp and H y lines. By calculat- 
ing a single line profile with an average SNR of about 1200 for 
each spectrum, we increased the original SNR of the individual 
spectra by a factor of six. The mean barycentric radial veloc- 
ity of the star is -10.8±0.2 km s _1 . We also derived a vsin ; 
value of 58+2 km s _1 from i) the first zero-point position of the 
Fourier transform for the mean LSD profile and ii) a nonlinear 
least-squares fit of an intrinsic profile computed for the stellar 
physical parameters of HD 50844 convolved with a rotationally 
broadened one. 

Before presenting the results of the frequency analysis, we 
clearly note that the shorter baseline of the spectroscopic obser- 
vations (27 d) implies a frequency resolution worse than that of 
the CoRoT photometric data. The power spectrum of the line 
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Table 3. First 20 frequencies identified in the amplitude spec- 
trum. 



Table 4. The frequencies identified in high-resolution spec- 
troscopy. 



Term 


Frequency 
[d- 1 ] 


Ampl. 
[mmag] 


Possible combination 
terms 




/i 


6.92528 


15.4539 


A 


,=2/l,/ 2 l 5 =3/! 




A 


11.21669 


6.7220 


fill = 


fl + A, /l23 = A - 


f 










fn - 2/2 




A 


11.25807 


4.0143 


flS9 = 


fi + A, A06 = A - 


f 


A 


12.84846 


3.3603 


fm = 


f + A> /791 = A - 


f 


A 


12.23831 


3.3271 


An = 


f + A, A95 = A - 


f 


A 


14.44689 


3.1767 




fig = f + A 




A 


13.27347 


2.6672 




As6 = f + A 




A 


13.35865 


2.3077 


/l49 = 


fi + fi> A 36 = fi- 


fi 


A 


11.75102 


1.4372 








Ao 


5.26653 


1.1277 




fe* = A + Ao 




fn 


14.46126 


1.0604 








J 1 2 


14.43209 


0.9731 








/l3 


9.95254 


0.9385 




fm = fi3 - A 




fu 


11.98380 


0.8268 




Al9 = fl + fu 




As 


6.55624 


0.8259 




/l41 = fl + As 




A, 


7.40446 


0.8099 








fn 


13.56880 


0.7768 




fmi = fn ~ A 




fia 


10.26062 


0.7332 




Am = fi + fm 




fl9 


5.78177 


0.7046 








Ao 


6.62914 


0.6435 




fm = f + Ao 





profile variations (LPVs) detected in the spectroscopic data is 
shown in the bottom panel of Fig.Q] Since the observations are 
single-site, the spectral window is more complicated because of 
the presence of relevant aliases at integer values of d -1 from the 
real peak. These facts result in uncertainties on the amplitudes 
and on the reliable detection of the spectroscopic frequencies. 
Nevertheless, the study of the LPVs is of paramount importance 
for corroborating and complementing the photometric results, as 
suggested by the different appearances of the photometric and 
spectroscopic power spectra (Fig.[TJ middle and bottom panels). 

The spectroscopic frequencies were searched for by 
analysing the LPVs by means o f the pixel-by-pixel technique 
dMantegazzal |2000t IZimal 12006b . Tableg]lists the terms detected 
in order of decreasi ng spectrosc opic amplitude (in continuum 
units as defined by IZimal 120061) . which is a measurement of 
the contribution of the given term to the line profile variabil- 
ity. We detected 27 terms with an SNR higher than 4.0 and 3 
terms with 3.1<SNR<3.5 (see Table 0). The three modes with 
3.1<SNR<3.5 were also considered in the LPV analysis since 
they are associated with well-defined C0R0T frequencies. The 
spectroscopic frequencies were associated with the closest fre- 
quency detected in all the photometric solutions. In the few cases 
in doubt (4 out of 30), the term with the largest amplitude was 
preferred. This assumption did not modify the subsequent {I, m) 
identification since the amplitude and the phase diagrams did 
not appreciably change with the frequency used to fold the spec- 
troscopic data. A different association with a photometric term 
would have modified the photometric amplitude, but this param- 
eter cannot be used in the spectroscopic mode identification. 

The mode identification was performed by fitting the am- 
plit ude and phase varia tions of each mode across the li ne pro- 
file dMantegazzal 120001) by using the software FAMIAS dZimal 
120081) . The identified {(,m) values are listed in Table [4] In the 
{(, m) couples, the negative m values indicate retrograde modes. 
Uncertainties are estimated to be ± 1 for the degree I and +2 for 



Term 


Frequency 


Spectr. 


SNR 


e,m 


EW RV 




[d- 1 ] 


Ampl. 






[kms -1 ] [kms~'] 


Ai 


14.43209 


0.27 


28.6 


4,2 


0.10 


Ai 


14.47568 


0.24 


22.2 


3,1 




A 


14.44689 


0.24 


23.0 


4.3 


0.009 


A 


13.27347 


0.18 


15.8 


3.2 


0.009 


A 


11.25807 


0.15 


13.1 


5.1 




A 


12.84846 


0.14 


19.5 


3,3 


0.008 0.52 


flA3 


15.54114 


0.14 


14.6 


9,2 




Ao 


15.22329 


0.13 


10.9 


12,10 




/is 


10.26062 


0.13 


12.5 


5,0 




fia 


13.85029 


0.13 


13.4 


2/, 




fm 


13.16891 


0.12 


13.0 


8,2 




Asi 


16.50451 


0.12 


9.2 


14,14 




A 


13.35865 


0.10 


8.9 


5,3 


0.11 


A 


6.92528 


0.11 


5.0 


0.0 


0.061 0.72 


A 


11.21669 


0.11 


9.5 


3.1 


0.022 0.39 


fa 


11.64327 


0,10 


9.3 


8.5 


0.08 


An 


15.12491 


0.10 


8.1 


7,4 




Ao 


5.26653 


0.09 


4.2 


2,-2 


0.08 


Ai 


15.80497 


0.08 


6.2 


4,3 




A 


11.75102 


0.08 


6.4 


11,7 


0.05 


fm 


19.71529 


0.07 


4.2 


12,8 




A* 


9.61191 


0.06 


4.1 


10,4 




fn 


5.67472 


0.06 


3.1 


4,-2 




An 


14.99570 


0.06 


6.7 


6,3 




Ah 


19.75025 


0.06 


4.1 


14,12 




A 3 


5.04299 


0.05 


3.5 


4,-2 


0.12 




5.49058 


0.05 


3.4 


5,-2 




A 4 


14.60061 


0.05 


4.1 


6,4 




A 


12.23831 








0.018 0.45 


fn 


14.46126 








0.20 


fn 


9.95254 








0.10 


Ao 


5.41862 








0.08 


fl9 


5.78177 








0.08 


Error 




0.02 






0.002 0.02 



the order m. They are larger in the cases of the modes /12 and, in 
particular, fn and fe, since they are very close in frequency. Also 
the nonsignificant detection of f\ \ is probably a problem related 
to the low resolution in the spectroscopic data. Mainly tesseral 
modes were detected, but also a couple of sectoral modes were 
found. This technique also returned a reliable estimate of the in- 
clination angle, i.e., ;=82+4 deg (Fig.[7]l. Therefore, HD 50844 is 
seen almost equator-on. The associated equatorial rotational ve- 
locity is not significantly different from the v sin i value derived 
above. This particular orientation explains why only one axisym- 
metric mode (i.e., /ig, Tabled was observed: these modes pro- 
duce marginal variations along our line of sight. 

An independent search for frequencies was performed on the 
timeseries defined by the th (equivalent width, EW) and 1 st (ra- 
dial velocity, (v)) order moments. Since these quantities are in- 
tegrated over the whole stellar disk, they are more sensitive to 
reveal low-degree modes. The amplitude s of the frequenc i es, de- 
tected in both moments by the Period04 dLenz & Bregerll2005l) 
method, are listed in Table [4] We recover in the EW and (v) 
variations (bottom of Table [4]i some of the terms with a large 
photometric amplitude that have been missed in the LPV analy- 
sis, because of their intrinsic low spectroscopic amplitude. Since 
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they were detected in the first moments, we can speculate that it 
probably concerns low-degree modes. 

The six frequencies detected in the EW variations corre- 
spond to six out of the seven first photometric frequencies. The 
term is probably missing because it is not well-resolved from 
fi and fs - 1 d _1 . The correlation between the terms found in 
the 1 st moment and those detected in the CoRoT photometric 
data is not as clear as in the case of the th moment. We note 
that the frequency 13.85029 d~'=2/i was detected in the pixel- 
to-pixel variations (Table @}. Does it really concern the first 
harmonic of f\ or is it an independent mode excited by reso- 
na nce? We tested the hypothesi s following the method described 
by iTelting & SchriiversT d 1997b . which states that a frequency 
is a first harmonic if the phase difference i/^i = 2if>o - i/ri be- 
tween the line-centre phases of the two frequencies is 1 .50 n rad. 
Because we measure 1//01 = 1-58 n + 0.19 rad, we conclude that 
13.85029 d~' is indeed the first harmonic of f\ . 

The main goal of the spectroscopic observations is to an- 
swer the crucial question about the huge number of frequencies 
detected in the CoRoT photometric data: is the large number 
of photometric terms related to the visibility not only of low- 
degree modes but also of high-degree modes? To investigate this 
we focus on the mode identification results, as derived from the 
high-resolution data, presented in Table|4] In Fig. |8]we plot the 
identified degree t in function of order m. This figure shows that 
the high-degree modes also tend to have a high order m. The 
highest detected degree is £ = 14. This observational fact tells 
us that we have at least 235 possible modes of pulsation for a 
given radial order n. A few different radial orders are sufficient 
for explaining the 1000 frequencies detected in the CoRoT pho- 
tometric timeseries of HD 50844. On the other hand, we know 
from CoRoT photometry that the frequency values higher than 
30 d -1 are observed after the 300 th rank of detection, i.e., with 
very low amplitudes. The observed spectroscopic frequencies 
have / <20 d -1 . We still miss the spectroscopic counterpart of 
the 20 < / < 50 d 1 region. Such a region could be filled by 
modes with high radial orders n, which are shifted toward higher 
frequencies than those of the low orders. 

The excitation of the almost totality of the modes allowed 
within a narrow radial order interval and its possible extension 
to a wide range of radial orders can explain the presence of 1000 
terms in quite a natural way. In either case, we can state that 
the huge number of photometric terms are plausible since the 
spectroscopic data are confirming that we are observing modes 
with high ((, m) values. In turn, this means that the cancellation 
effects for modes up to t ~ 14 are not working and that these 
modes are able to produce a net flux measurable at the 10~ 5 mag 
level. 



4.2. High-resolution spectroscopy: abundance analysis 

In addition to the LPV analysis, the FEROS spectra were used 
to determine the abundances of the elements. The atmospheric 
mode ls were compu ted with the line-blanketed LTE ATLAS9 
code dKuruczL Il993l) . which handles line opacity with the opac- 
ity distribution function. The synthetic spectra were computed 
with the SYNTHE code dKuraczl [19931) . The stellar-line identi- 
fication and the abundance analysis were perfo rmed on the ba- 
sis of the VALD lines list dKupka et all l2000b . The mean LTE 
abundances log e(El) (by convention, log e[H] = 12), along with 

the uncertainties determined as <x 



100 



t + "Toe b are reported 



in Table The atmospheric parameters T e ff = 7400 K, log g = 
3.6 and £ = 2 km s _1 were adopted. To determine the uncer- 




40 60 80 

Inclination [deg] 

Fig. 7. Behaviour of the x~ values, obtained by fitting all the 30 spec- 
troscopic modes to the LSD profiles, each time keeping the same fixed 
value for the inclination angle. 



15 



10 



a 5 



-5 



n — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — r 

A 



A 



A 

A 



A 



A 

AAA A 
▲ ▲ ▲ A 

A A A A 

A 

A A 
A ▲ A 



_i 1 1 1 1 1 1 1 1 1 i_ 



10 



15 



L 



Fig. 8. Identification of the (€, m) couples from high-resolution spec- 
troscopy. Open triangles: one term with that (f, m) identification. Filled 
triangles: two terms with that (£, m) identification. 



tainties connected with effective temperature and surface grav- 
ity, we calculated the abundances for atmospheric models with 
r eff = +200K and logg = + 0.2. 

As a general result, HD 50844 is a slightly metal-deficient 
star. However, Table [5] shows how the abundances of elements 
as C, N, and S in the HD 50844 atmosphere are very similar to 
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Fig. 9. Element abundances of HD 50844 derived f rom FEROS spectra 
compared to the solar values dGrevesse~et "aTl l2007l) . 



Table 5. Mean LTE abundances log e(El) along with the uncer- 
taintie s, cr, with solar abundances determined bv lGrevesse et al.l 
(I2007I) . 



Element 


Z 


loge(El) 


<r 


Sun 


C 


6 


8.45 


0.130 


8.39 


N 


8 


8.72 


0.102 


8.66 


Na 


11 


5.75 


0.097 


6.17 


Mg 


12 


7.02 


0.113 


7.53 


Si 


14 


6.95 


0.386 


7.51 


S 


16 


7.25 


0.054 


7.14 


Ca 


20 


5.91 


0.112 


6.31 


Sc 


21 


2.56 


0.103 


3.17 


Ti 


22 


4.40 


0.159 


4.90 


V 


23 


3.38 


0.146 


4.00 


Cr 


24 


5.10 


0.140 


5.64 


Mn 


25 


4.68 


0.157 


5.39 


Fe 


26 


6.71 


0.126 


7.45 


Co 


27 


4.89 


0.098 


4.92 


Ni 


28 


5.53 


0.151 


6.23 


Cu 


29 


3.52 


0.018 


4.21 


Zn 


30 


3.93 


0.079 


4.60 


Y 


39 


1.67 


0.099 


2.21 


Zr 


40 


2.37 


0.458 


2.58 


Ba 


56 


1.91 


0.170 


2.17 


La 


57 


0.483 


0.126 


1.13 


Ce 


58 


0.292 


0.436 


1.70 



nal uvby dataset consists of 1496 datapoints, collected in 66 
nights for a total survey of 317.4 hours. Unfortunately, the space 
and ground-based observing windows are not simultaneous. The 
spectral window of the ground -based data is characterized by 
several alias peaks (at integer values of d -1 and ys~' of the cen- 
tral peak) and the true peaks could only be identified knowing 
the CoRoT solution. We considered the first 20 CoRoT frequen- 
cies and refined their values by performing a least-squares fit of 
the data in the v colour, i.e., the one with the largest amplitude. 
This step was necessary since the time baseline of the uvby data 
is much longer than that of the CoRoT data. The first 13 fre- 
quencies resulted in a significant v-amplitude. The v-amplitude 
of the f\ term is 22.2+0.3 mmag, while those of the others are 
much lower (e.g., 8.0 mmag for/2 and 1.70 mmag for /13). 

We calculated the amplitude ratios an d the phase sh ifts by 
using the y colour as the reference system (iGarridol l2000h . They 
were compared with the theoretical predicti ons computed us- 
ing th e non-adiabatic pulsational code GraCo dMova & GarridoL 
120081) . The immediate result is that the predominant term 
/i=6.92 d _1 is characterized by positive phase shifts (i.e., 4> u , v ,b- 
<p Y > 0.0). This puts a very tight constraint on the mode identi- 
fication, since only the radial (^=0) mode shows positive phase 
shifts in the (pi, - <p y case (Fig.fTOt. On the other hand, the com- 
parison between the observed and the calculated values does not 
supply conclusive mode identifications in the cases of the other 
terms. 

To further investigate the f\ mode identification, we note 
that the harmonics 2fi=fzg (amplitude 0.56 mmag) and 3/i=/2is 
(0.12 mmag) have been detected in the CoRoT data. The ob- 
served values differ from the exact multiples of f\ by 3 x 10~ 4 
and 1 x 10~ 3 d _1 , respectively. Both the 021=4.04 rad and the 
031=2.11 rad Fourier parameters are very similar to th ose ob- 
serve d in high amplitude 5 Set stars (Figs. 4 and 5 in iPorettil 
|2001|) . known to pulsate in the fundamental (rarely in the first 
overtone) radial mode. Identification of f\ as a radial mode is 
also consistent with being the highest amplitude term both in ra- 
dial velocity and in equivalent width. Although a nonradial mode 
could display such large variations, such variations are usually 
connected with the large displacements typical of radial modes. 
The mode identification of the f\ term from the spectroscopic 
data was not trivial since we had to take the harmonic 2/j and 
the EW variations into account. The FAMIAS method returns 
two possibile (t,m) couples, i.e., (0,0) and (2,0). When consid- 
ering both the spectroscopic and the photometric identifications, 
we can consider the f\ term as the fundamental radial mode. 



those of the Sun, while other elements are underabundant (see 
also Fig. Pi. Suc h an abundance pattern is typical of A Boo stars 
JPaunzenl |2004|) . Wi th respect to the claim ed spectroscopic bi- 
narity of these stars dGerbaldi et all 120031) . we did not find any 
evidence of a second component in the FEROS spectra, and we 
did not detect a long-term radial velocity drift that could be as- 
cribed to an orbital motion. 



4.3. Multicolour photometry 

The importance of colour in formation to es timate the degree £ 
of the modes is well-known (Garrido, 2000'); therefore, we ob- 
served HD 50844 in Stromgren uvby photometry to accompany 
the CoRoT white-light photometry. Several campaigns were car- 
ried out at S. Pedro Martir and Sierra Nevada observatories, 
using twin Danish photometers. The multicolour observations 
span 743 d. The campaigns started when the IR01 pointing was 
decided and continued after the CoRoT observations. The fi- 



5. The seismic inferences 

The discussion of the number of excited frequencies in 6 Set 
stars has followed a particular pattern. At the beginning of the 
observa t ional efforts to disc over as many modes as possible 
(IPorettil l2000t iBregerl 120001) , the number of detected frequen- 
cies was less than expected and a kind of sele ction effect in excit- 
ing modes/dampin g amplitudes was c l aimed ( Guzik et aQfeOQi 
and, more recentlv. lBreger et al.l2008l) . lGarrido & Porettil (1200- 
suggest that the problem was mainly of an SNR nature. They 
stressed how the number of detected terms in the science case 
of FG Vir follows an exponential behaviour as a function of the 
observational technique used, namely from single-site observa- 
tions in a few nights to multisite observations in several years. In 
a certain sense, they correctly predicted that a long-term, high- 
precision space mission as CoRoT would allow access to the 
distribution of the amplitudes down to 10~ 5 - 10~ 6 mag, allow- 
ing the harvest of a multitude of modes. That the distribution of 
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Fig. 10. Comparison between theoretical predictions and multicolour 
photometric observations for /j , the mode with the largest photometric 
amplitude. Amplitude ratios vs. phase differences (in degrees) in b and 
y colours are shown. Red large square is the observed errorbox. Points 
represent theoretical predictions for spherical degree ^=0 (black filled 
square), C=l (blue star), and (=2 (green open square). 
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the frequencies is quite uniform in the 0-20 cT 1 region (Fig. |4j 
implies that modes are excited at lower frequencies than that of 
the fundamental radial mode, i.e., / <6.9 d _1 . The peaks with 
f < 3d 1 can be enhanced by an incomplete removal of the in- 
strumental drift. The spectrum of HD 292790 (Fig. [3] top panel) 
shows a noise level below 0.05 mmag in this region and there- 
fore this contamination should become important only after the 
600 th frequency (see Table [TJ. We can argue from the CoRoT 
timeseries on HD 50844 that gravity and mixed modes are quite 
common in 6 Set stars. 

The frequency set obtained from the CoRoT timeseries of 
HD 50844 allows us to investigate the problem of the regu- 
lar spacings in the frequency spectrum of 5 Set stars. There 
have been several claims that such spacings will be modified by 
nonlinear effects and m i xed modes, especially in evolved stars 
dChristensen-Dalsgaardl 120001: iGoupil & Dziembowskil 120001: 
iGuzik et al.Ll2000t) . This hint is confirmed by the frequency spec- 
trum of HD 50844, which does not supply any reliable indica- 
tion of the regular s pacings detected in solar-like pulsators (e.g., 
iMichel et aUl2008h . 

As a first attempt to model the scientific case of HD 50844, 
equilibrium models representative of the star have bee n calcu- 
lated with the evolutionary code cesam dMorell 11997b . taking 
first-order effects of rotation into account. Such models (the so- 
called pseudo-rotating models) include the spherically averaged 
contribution of the centrifugal acceleration by means of an ef- 
fective gravity. The nonspherical components of the centrifugal 
acceleration (not considered in the equilibrium models) are in- 
cluded in the adiabatic oscilla tion computa t ions b y means of 
a linear perturbation analysis dSuarez et a 120061) . Figure QT| 
shows the models within the uncertainty box in T e ff and log g. 
The models, evolved considering the local conservation of the 
angular momentum, are able to return a surface rotational veloc- 
ity around 60 km s _1 (i.e., the value obtained from vsim taking 
; = 80 deg as suggested by the analyses of the the line profile 
variations) starting from a value of 75 km s _I on the ZAMS. We 
can also compute the fundamental radial mod e of our models us- 
ing the code GraCo dMova & Garridoll2008h : the darker box in 
Fig. [TTI marks the cases matching the observed value f\. When 
considering the slightly metal deficient content of HD 50844, 



Fig. 11. HR diagram containing several evolutionary tracks of 
representative models of the star. The shaded area represents the 
observational uncertainties in effective temperature and gravity. 
The models in the darker area have a fundamental radial mode 
equal to f\ =6.925 d _1 . The filled star represents the observa- 
tional average location of HD 50844 in the diagram. Continuous 
lines represent tracks of evolved models with a solar metal- 
lic content, whereas dotted lines represent evolutionary tracks 
of submetallic models; the latter are more representative of the 
HD 50844 case. 



we can argue that the models suggest a mass between 1 .65 and 
1.80 M Q . 

6. Conclusions 

The exploitation of the CoRoT photometric timeseries of 
HD 50844 resulted in a very complex and intriguing task, 
also giving a completely new picture of the pulsational con- 
tent of 5 Set stars. We demonstrated that the light curve of 
HD 50844 can be explained by the presence of hundreds of ex- 
cited terms. Classical checks such as using different software 
packages, subdividing the timeseries into different subsets, care- 
fully inspecting the residuals, and comparing with other sim- 
ilar targets observed by CoRoT, together confirm this issue. 
High-resolution spectra support the excitation of modes hav- 
ing a very high degree ( (up to £=14), thus providing an ob- 
servational explanation for the richness of the frequency spec- 
trum. Such a large number of modes of different £ simultane- 
ously excited is among plausible features sugge sted by theoret- 
ical works (e.g.. lBalona & Dziembowskil, 1 19991) . An immediate 
conclusion is that the cancellation effects are not sufficient in 
removing t he variations of the flux integ r ated o ver the whole 
stellar disk; iDaszvhska-Daszkiewicz et all (120061) correctly pre- 
dicted this extension to high degrees of the modes extracted from 
high-precision photometric timeseries. We should keep in mind 
that the CoRoT performances and the continuous monitoring al- 
low detection of amplitudes of about 10~ 5 mag without signifi- 
cant aliasing effects, which is a totally new perspective. 
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The case of HD 50844 seems to match the theoretical picture 
featuring a very large number of excited modes with amplitudes 
limited by a saturation mechanism of the k mechanism. This con- 
firms that the observation of 6 Set stars with CoRoT will allow 
us to address the longstanding problem of amplitude distribution 
and the potential existence of mode selection processes in the 5 
Set instability strip. We note that HD 50844 probably does not 
constitute a "standard" scientific case of a 6 Set star, because it 
is located on the TAMS, probably belongs to the class of A Boo 
stars (i.e., showing atmospheric particularities), and it is seen al- 
most equator-on. Nevertheless, seismic models that were able to 
reproduce the position of HD 50844 in the HR diagram and its 
measured rotational velocity were calculated. 

HD 50844 has for the first time disclosed the intrinsic com- 
plexity of frequency spectra of intermediate-mass stars in the 
lower part of the instability strip. It is quite evident that, after 
these first CoRoT observations, we must look at 6 Set stars in a 
completely new way. We have to realize that ground-based ob- 
servations were only able to observe the tip of the iceberg, while 
the larger part of the modes have remained undetected. 
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